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Abstract 

Within the theoretical framework of divergence- type theories (DTTs), we set up a consistent 
nonhnear hydrodynamical description of a conformal fluid in flat space-time. DTTs go beyond 
second-order (in velocity gradients) theories, and are closed in the sense that they do not rely 
on adiabatic expansions. We show that the stress-energy tensor constructed from second-order 
conformal invariants is obtained from the DTT by a consistent adiabatic expansion. The DTT 
satisfies the Second Law, and is causal in a set of fluid states near equilibrium. Finally, we 
compare, analytically and numerically, the equations of motion of the DTT and its truncation to 
second-order terms for the case of boost invariant flow. Our numerical results indicate that the 
relaxation towards ideal hydrodynamics is significantly faster in the DTT than in the second-order 
theory. Not relying on a gradient expansion, our findings may be useful in the study of early-time 
dynamics and in the evolution of shock- waves in heavy- ion collisions. 
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I. INTRODUCTION AND MOTIVATION 



There is currently a great interest in relativistic dissipative hydrodynamics, mainly due 
to its application to the description of the hot dense QCD matter created in the Relativistic 
Heavy Ion Collider (RHIC) experiments 
The application of the AdS/CFT correspondence 
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study strongly coupled conformal plasmas unaccessible to kinetic theory, has also fueled 
considerable interest in the formal aspects of this formalism. The study of the hydrodynamic 
regime of conformal field theories is important since QCD is approximately conformal at high 
temperatures [25 1. 



The need for dissipative corrections in modeling heavy-ion collisions is (at least) two- 
fold. First, quantum uncertaint y p revents the existence of a perfect fluid. Since first- 
order relativistic hydrodynamics [261 . 0] is known to have strong drawbacks, among them 
lack of stable solutions and acausal propagation of perturbations, one should really go to 
second-order theories (see however Refs. [28|, l29|]). Second, the description of heavy- ion 



collisions in terms of perfect hydrodjTiamics works well in almost central Au+Au collisions 
near midrapidity, but gradually breaks down in non-central collisions and at forward rapidity 
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15|. 



The complete second-order stress-energy tensor of a strongly coupled conformal fluid in 
d = 4 space-time was given independently by Baier et al 20|], and by Bhattacharyya et al 21|. 
Recently, Loganayagam [sO] developed a very useful Weyl-covariant formalism and proposed 
a local entropy current consistent with the second-order T„,^ derived previously and with the 



3l|). One of the most important 



Second Law (see also related work of Romatschke, Ref. 

results of Refs. 20, 1211, 0] is that these works show that the hydrodynamic description of a 



33| formalism 



conformal fluid does not belong to the conventional Israel- Stewart 
(see also Refs. for a discussion of this issue). This is because the conventional (or 

entropy- wise) IS theory, not being a controlled gradient expansion 0, 14|, cannot account 



'or shear-shear coupling, which is p resent in the stress-energy tensor of the conformal fluid 
3I, y, 0, 31] (see Refs 



• [14, ll5|, [16(1 for interesting discussions on this and related issues in 



the context of dissipative fluid dynamics as derived from kinetic theory). 

The main purpose of this paper is to set up a consistent hydrodynamical description of a 
conformal field theory within the theoretical framework of divergence-type theories (DTTs) 



3 



34 l35l . l36l . l37l . l38l . l39l . |40|]. Our goal is to go beyond second order theories (in velocity 
gradients) by formulating a theory in closed form, that is, without reliance on adiabatic 
expansions. We will not tackle the full problem of a conformal fluid in curved space-time 
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3l| . but limit ourselves to Minkowski space-time. Another aim of this work 



is to analyze the causality properties of the DTT developed here, and to compare our 
results for the stress-energy tensor with those obtained from the derivative expansions of 
Ref. 3, [21I, We note that we do not calculate transport coefficients in this paper, but 
assume they are known either via kinetic theory or the AdS/CFT correspondence. 

DTTs are interesting alternatives to the IS formalism (although they may be physically 
equivalent in certain cases) because the conditions for hyperbolicity and causality of the full 
nonlinear evolution can be stated in very simple terms. As clearly shown by Liu, MuUer 
and Ruggeri 36|], DTTs are often more general and flexible than the IS theory, allowing a 



systematic derivation of nonlinear terms in constitutive equations, which are not captured 



by entropy-wise IS theory (see Refs. 
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16, 



20 



23(1). Besides, they have the extra 



advantage that, being the equations of motion of divergence type, discontinuous solutions 
(shocks) can be given mathematical meaning, which is relevant to the phenomenon of conical 



flow in heavy-ion collisions 
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24l |. Moreover, and this is an important point for what 



follows, the symmetries of the theory can be coded directly on the generating function of 
the DTT. 

The main results we arrive at are that: (i) for the case in which the second-order transport 
coefficients A2 and A3 vanish, the second-order (in velocity gradients) stress-energy tensor 



constructed from conformal invariants 
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30|, can be consistently derived via an adia- 



batic expansion from the DTT we set up; and (ii) the DTT and its adiabatic expansion are 
causal for states near equilibrium, and satisfy the Second Law. We also obtain, as a simple 
illustration, the hydrodynamic equations of the DTT for the case of boost invariant flow, 
and compare them to those of the second-order theory. Our numerical results show that the 
DTT approaches the ideal fluid behaviour faster than the second-order theory. 

We believe that the DTT presented here may be useful in the study of two aspects 
of heavy-ion collisions, both of which seem to require theories going beyond second-order 
velocity gradients. First, early-time dynamics, where velocity gradients are not small and 
for which even IS formalism shows unphysical behaviour such as reheating |l^ (see also 
Chesler and Yaffe 2J], who study the creation and evolution of a boost invariant anisotropic 
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plasma directly from the gravity side, and El, Xu and Greiner 111, w ho develop a novel third 



order theory). Second, the evolution of initial state fluctuations [38j, for which higher order 
terms are crucial (see particularly the work of Lublinsky and Shuryak 2^, who developed 
a linearized hydrodynamical theory that includes, in principle, all-order velocity gradients). 
In this respect (that of "resumming" higher order velocity gradients), the developments of 



Ref. [23j are related to ours. We note that the DTT developed here containts, in addition, 
all quadratic terms in velocity. 

The paper is organized as follows. In Section [Tll we briefly review divergence-type theo- 
ries. In Section HIT] we first review some basic properties of the hydrodynamics of conformal 
field theories, and then set up the divergence-type theory of a conformal fluid. We also prove 
that the DTT satisfies the Second Law exactly, and obtain the hydrodynamic equations. In 
Section [IV] we investigate the causahty properties of the DTT, for fluid states near equilib- 
rium. In Section |Vl we show that the dissipative part of the stress-energy tensor obtained 
from second-order (in velocity gradients) conformal invariants can be obtained from a con- 
sistent adiabatic expansion of the DTT. We note that the DTT cannot reproduce terms 
containing the vorticity tensor. In Section |Vl] we compare, analytically and numerically, the 
hydrodynamic equations of the DTT and of the second-order theory for the case of Bjorken 
flow. The paper closes up with a brief summary of results. 

II. DIVERGENCE-TYPE THEORIES 

In this section we give a brief summary of divergence-type theories (DTTs). Detailed 



discussions can be found in Refs. 



34 



38U39|,|40( (see also Ref. |36|). 



According to Geroch and Lindblom [34], the hydrodynamical description of a nonequilib- 
rium state requires, besides the particle current Na and the stress-energy tensor T^^, a new 
third order tensor A^^p obeying an equation of motion of divergence type. The dynamical 
equations are the conservation laws of A*"^ and T^j^, together with an equation describing the 
dissipative part: 

A^'^P^ = FP (1) 

where A and / are algebraic local functions of and T and symmetric in the indices (i/, p). 
A semicolon stands for a covariant derivative. The entropy current is extended to 

Sf" = <^f' - I3^T'"' - aA^^ - A^'^P^^p (2) 



where (5^ = Uy/T is the temperature vector, a = /i/T is the affinity, $^ is the thermodynamic 
potential and C,up is symmetric, traceless and vanish in equihbrium. Note that, in equihbrium, 
13'^ is Kilhng and a is constant. 

We now require that the entropy and the thermodynamical potential be algebraic func- 
tions of {a, P^^^^y). If the entropy production is to be nonnegative, then 

/}$/^ ^(^^^ 

V = ^"; ^ = ^ = A^'^p (3) 

Thus, as a consequence of the equations of motion, the entropy production rate is 

s^, = -r%p . (4) 

Since the stress-energy tensor is symmetric, we must also have 

where xi'^y f^^ly ^fj.u) is the so-called generating function of the theory. This means that 
every DTT is completely determined once x ^i-nd / are specified as algebraic functions of 
a, P^, ^^u- The theory thus constructed satisfies the principles of relativity and entropy, and 
fully exploits the latter jscl . 

Introducing the symbol to denote the set (a,/?^,^^^), the set {N^' ,T'"' , A^"'P) and 
Ib the set (0, 0, 1^,^), the theory is summed up in the equations 

The equations of motion can also be written as 

M'bcC = (7) 

where 

M^^ = M^^ = a^^VaC^aC^ . (8) 

The system of equations ([7]) is automatically symmetric since matrix M'^^j is symmetric in 
the indices A and B due to the fact that partial derivatives conmute. Causality is therefore 
ensured if the quadratic form M'^^uj^ is negative definite for all future directed timelike 
vectors cj^, or equivalently if = Mp^SC^SC'" is timelike and future oriented for any 



displacement from an equilibrium state 3J, |35|. Note that causality depends on the 

BC 



form of Mo^ and not of Ib 



III. CONFORMAL HYDRODYNAMICS AS A DIVERGENCE-TYPE THEORY 



We consider the hydrodynamic regime of a conformal quantum field theory in d = A flat 



space-time (for a general discussion, not limited to conformal fields, see Ref. [?, 4l|)- In such 
a theory, the classical action evaluated on the classical equations of motion is invariant under 
a Weyl transformation g^,^ — > e~'^^^^^^g^,y, where lu is a function of space-time coordinates 
x'^. The classical stress-energy tensor of such a theory is necessarily traceless, while the 
quantum one presents Weyl anomaly. As shown in Ref. 20|, in even dimensions d the 
number of derivatives appearing in the Weyl anomaly is precisely d, which means that 
second-order hydrodynamics in d = 4 dimensions is Weyl invariant. It is straightforward to 
show that, for a conformal theory. 



(9) 



under a Weyl transformation. Therefore, for a conformal fluid the conservation law of energy- 
momentum T.^'^ = 0, where a semi-colon denotes covariant differentiation, is automatically 
Weyl covariant (see, for example, Refs. j^, 0, 42]). 



The tracelessness condition T^^ = imposes p = {d— l)p and C = 0, where p is the energy 
density in the local frame, p is the thermodynamic pressure, and ( is the bulk viscosity. 
The transformation rule for T^^" implies p e'^'^p, the four-velocity ■u'* e^u^ and the 
temperature T e'^T, which means that the temperature vector jS'^ = u^/T has conformal 



weight equal to zero. In addition, Eq. (Ml i mp 



A a constant (see, for instance, Refs. 



21 



ies that the shear viscosity r] = AT'^ ^, with 
3l|). In equilibrium P{fj.;u) = 0, which 



means that is a Killing vector (parenthesis around indices stand for symmetrization) . 
The four-velocity is normalized as u^u'^ = —1, and we use the signature (— , -|-, -|-, -|-). We 
will make use of these properties in what follows. 

A divergence-type theory is completely specified by its generating function x(a, /J'^, ^'^'^) 
and source tensor Ip^{a, j3^,C,^'^), where {a, are the fugacity, the temperature vector 

and a symmetric and traceless tensor that vanish in equilibrium, respectively. Our start- 
ing point will therefore be the specification of x ^-^^d / as algebraic functions of the set 
(a, ^'^''). We will deal with the source tensor later on, so for the moment we focus on x- 

Since we want to construct a quadratic DTT (in deviations from equilibrium), we will 
consider terms which are at most quadratic in the nonequilibrium tensor C,^'^. For simplicity. 



7 



we will restrict ourselves to a conformal theory with no conserved charges, which means 



a = 0. In this case, it is convenient to employ the energy or Landau-Lifshitz frame 
The generating function which satisfies these requirements can be written as 

X = X{0) + X{i) + X{2) 

3 



(10) 



with 



per 

1=1 



(1) 3 



(11) 



This is the most general local scalar constructed from T, g^u, and ^^i,, which is quadratic 
in the latter. In Eq. f|TT]) . A^'^ is the spatial projector orthogonal to u^. From the trans- 
formation rule for g'^'^ and u'^ we immediately obtain A'^'^ ^2uj^ij.u^ rjj^^ tensors S^^'"^ 
produce the most general decomposition of a symmetric and traceless tensor around a time- 
like direction u'^. 

As already mentioned, conformal invariance requires T^^'^ —>■ ^i.'^+'^)^T>^'' and = 0. In a 
DTT, this means 



9 X , ,id+2)u d X 



6^'^+^^'^ ^ (12) 



and 



Note that 



dp. 

= T^[/\^^j3^ + 13" A^^] and (14) 



3,4 

In the following, we will describe how these conditions determine the scalar functions xo,i 
and appearing in Eq. ffTOj) . It is clear that the conditions given in Eqs. f fT2|) and f|T3|) 
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must be satisfied separately by the zeroth, first and second order terms of the expansion of 
X- This is because truncating the expansion at zeroth, first and second order does not break 
the conformal invariance of the resuhing hydrodynamic theory. In other words, the zeroth, 
first and second order terms of T^^ are independent of each other. 

A. Perfect fluid 

At zeroth-order, condition (fT3l) on the stress-energy tensor 

T^- = pg^^- + u^u'Ip + p] (15) 

implies p = {d — l)p. In a DTT we have (recall that p and p are equilibrium quantities and 
thus completely determined by xo) 

To'' = TV'^ + T W (3^ + T^) , (16) 



which implies 



p = T — and 



Therefore, we find that in a conformal DTT xo must satisfy 

d'^Xo ^ jd - 3) rfxo 
c^T2 ~ T dT ' 

The solution to Eq. (ITS]) is 

Xo = aT'^-^ + a' (19) 

where a and a' are constants. Note that a' is irrelevant since it does not change T^^, so we 
set it to zero. 

From the transformation rule for T we immediately obtain xo ~^ e^'^~'^^^Xo- Taking into 
account that /5° is Weyl invariant while 



dxo , ^d'^Xo 



(17) 



(18) 



- ^ = e-^-/3« (20) 



we have 



' -8(3,8(3.^' 

_ Jd+2)ulrplJ.U 

— e Jq 



(21) 
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as it should. 

From the above it is clear that, in order to obtain the correct conformal weight for the 
complete second-order stress-energy tensor (Eq. fll2p ). the generating function given in Eq. 
fITU]) must transform like Xo under a Weyl transformation. That is, we require that 

X^e^'-'^-X . (22) 

As already noted, this implies that the first and second order terms ^ud X{2) have 
conformal weight equal to {d — 2). 

B. Linear DTT 

We will now determine XiW ^-ud from it the stress-energy tensor at first order in the 
nonequilibrium tensor ^a/s- 

rpt^y ^ rj.^^ ^ (23) 

with 

^^lv _ 9 X{i) 



It is convenient to rewrite the first order term in x as 

= Xi^p.u'W' = MpaPPp" (25) 

with Xi = T'^Xi- The tracelessness condition g^^r^" = implies the following differential 
equation for xi- 

whose solution is 

Xi{T) = h + cT^ (27) 

with h and c constant. 

The issue of how to choose the integration constants h and c is not trivial. The criterium 



of a bounded solution is not compelling enough |44l| to support the choice of constant c = 0. 
However, we would like to point out that the requirement that the generating function must 
have a definite transformation law under conformal transformations means that the two 
constants h and c cannot be both different from zero at the same time. Therefore there are 
two families of conformal divergence type theories: one with 6 = 0, c 7^ and the other 

10 



with b ^ and c = 0. Moreover, the physical content of both familes is the same. Indeed, 
if we substitute b by cT^, but after computing the energy- momentum tensor we replace 



fJ,U 



by T~^^^i, we obtain once again Eq. (HT]) . Of course, to reach this conclusion we use that 
the physical nonequilibrium tensor ^fj_y is transverse with respect to the four-temperature 

(this fact will be shown in what follows), so extra terms in the energy-momentum tensor 
vanish identically. Given the physical equivalence of both classes of theories, we have chosen 
to investigate only the 6 7^ matter of simplicity. 

As it can be seen from Eqs. fl22l) and (!25l) and the fact that (3^ is Weyl invariant, setting 
c = implies that the conformal weight of ^fj,u is equal to d — 2. Therefore, 

^au ~^ e^*^ "^^^au and 

(28) 

From the generating functional determined above, the tensor of fluxes 

- fife 



becomes 



where we have used that 



A-^^T = 6(/"/3-> + ^■>'5/?") - V^"^ , (30) 

2 



1^ = ^ (^p ^ J + ^:^;) - \9,j''' (31) 

which follows since S,^u is symmetric and traceless. 
The divergence of the tensor of fluxes is 

^5a7^ _ ^^(a;7) _ 5«7 . (32) 

The first order divergence-type theory is summed up by Eq. fl32l) together with 

rf = bi^"^ (33) 

and 

A'^% = /°^ . (34) 

The system of equations 03211341] must lead to the same r'^" of Eckart's theory (for a 
conformal fluid), which can be written as 



34 



36 



38|: 



- -vs^r («(p;.) - ^^(p) = -^T^rr^p;.) (35) 
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where the last equahty follows from the transversality of S^^-^ . Note that, since 

(36) 

•^(l) ^ "^(1) ' 

we must have, from Eq. (!35|) . that r] e^'^t], and therefore ri{T) = const. T^. 

In order for the first order stress-energy tensor obtained from the DTT to coincide with 
that of Eckart's theory we must provide a linear relationship between the source tensor I""^ 
and the nonequilibrium tensor 

= _/)«7P^^^^ _ (37) 
Using Eqs. (EUSD and Eq. ([351) we obtain 

where we have used that 

g,,S^;'{'' = A^(^A;) - ^A^^AP- = (39) 

since 

= ^ ^al3 ^ ^/Sa A"^ = 3 . (40) 

Since the conformal weights of ^po- and are both equal to 2 (see Eqs. fl32|) - flMll ). it 
is seen from Eq. (1371) that the tensor D'^'^^'^ must be Weyl invariant. From Eq. (138!) this 
implies r/ oc T^, as before. Note also that the requirement t^^ = r("^ automatically implies 

ipupa 



Pp.^^'^ = 0, since 5*^^^^ is transverse. We have 

where cTq,^ = Sii^a-ypuU^'^'^^ ■ The physical meaning of ^'^^ being transverse is that the bulk 
viscosity and the heat flux, which are both proportional to /5°^a^ [38] , vanish. The vanishing 
of the heat flux is expected since the chemical potential is zero, whereas the bulk viscosity is 
zero since the theory is conformal. In the next section we will show that the transversality 
of ^^'^ holds in the quadratic theory as well. 
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C. Quadratic DTT 



We now go over to the quadratic stress-energy tensor given by 

T'^^ = T^'' + rr + rr (42) 

with 

From the conformal weights of given in Eq. ( j28l) . and of S'^^'"^ (see Eq. ( l36l) ) it is 
seen that, to obtain the correct conformal weight for T^*^, we must have X2 = CiT~^, where 
Ci are constants to be determined. This ensures that X{2) in Eq- ( flOl^ has conformal weight 
equal to (c? — 2), which means r^'^ has conformal weight equal to {d + 2). 

The tracelessness condition, Eq. f|T3l) . will determine relations among the coefficients q. 
The quadratic contribution to the stress-energy tensor can be written as 

rr = r'^^X:Q^(7'^e«.U (44) 
1=1 

where we have defined the operator 

_d . d 
d 



In this notation, the trace of r^'^ becomes 

3 

9,uTr = T^,Y.^.S-^'"^^'^,^,. (46) 
1=1 

with 

Computing the derivatives and equating the coefficients of the Lorentz invariants 

A^^'A^'^e.^epa, P'^P^A'^^^^Cp. and P'^P^ P^P'^^a.^p. (48) 

to zero {g^vT2^ = 0), we find a linear system of equations for the three unknowns q. It turns 
out that the equations involving C2 and C3 (which come from the invariants /5^A^^ and /?^^^) 
are inconsistent with each other, the only way out being imposing that 

= . (49) 
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That is, the tracelessness of forces the transversahty of the nonequihbrium tensor ^ap- 
This means that the heat flux and the bulk viscosity remain zero at second order, which 
is a satisfying result. With this additional requirement on the nonequihbrium tensor, C2 is 
left unspecified, while the remaining equation (coming from the invariant A"'^A''''^^„^^po- and 
relating Ci and C3) reads 

ci = -^C3 . (50) 

Therefore, we have found that the quadratic part of generating function of a conformal 
fluid can be written as 



X(2) 



T- 



I qa'ypu 



(51) 



From Eq. (!5T!) we can calculate the tensor of fluxes 



2) 



(52) 



a7 



where A^^"' is the first order term given by Eqs. (!29l) and ( |30l) . We can rewrite A2 as 



A 



per ) 



(53) 



with 



+ 2c2 A'^^ ( AP'^P^ + A^^ 



+ 2c2 A^^ |^A""/3^ + /5"A"^^ 



The divergence of ^12"^ is 



rin^Oi-ypa 

= %;^P^;si 



^Scx'jpa 

d^pa 



(54) 



(55) 
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We get, after some algebra, 



+ Kr^'"-'^ i,„P^,s + G'^''''ip.;8 (56) 



with 



In order to have a complete theory at second order in deviations from equilibrium, we 
must find a suitable source tensor I2 quadratic in We will find the constraints imposed 
on I2 by requiring that the Second Law holds, and by the fact that C,^"^ is traceless and 
transverse, and find an explicit expression for 12- Guided by linear results, we will consider 
that I2 has the form 

ir = J''^'"""'^pa^,. , (58) 

where J = J(/3^ A^"). 

In a DTT, the entropy production is simply 5'.'^ = — /"'''^q,^. We have 



CA* _ CM 



1 



2 



(59) 



S/^li is the entropy production of the linear DTT (in ^), and is clearly positive definite. The 
problem comes from 5'.'^|2, which has no definite sign. In order for the Second Law to hold 
for arbitrary ^, we must require that 

ir^a, = . (60) 

We will now find the explicit form of J^tpo-/^^^ (ggg Eq_ (j^Sj)). Since is a local function 
of jS^ and A^"^ , ^^'^ is traceless and transverse, and Eq. fl60|) must hold, we see that J"')'^'^^'^ 
can only have two terms, one proportional to (3'^I3"' and the other to A"'*'. In addition, the 
only non- vanishing scalar we can form out of ^pa^^v is 

Sllfipai,. = i'-'ip. . (61) 
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Therefore, we have 

ja^pai^u ^ ^y.^ (T)/5°/5^ + /2(T) A"-^^ Sl^^^" (62) 

where the fi{T) are functions of temperature we must determine. It can be checked from 
Eq. ([56]) that 

/?„/?^Af7 = , (63) 

which means that /i = 0. In order to find f2{T), we will consider a power law dependence 
and use the fact that I2 has conformal weight equal to 2. Recalling that A°''^ and S^^^'^'^ 
have conformal weights 2 and 4, respectively, we immediately obtain /2 = gT~^, where g is 
a constant. 

So, the final expression for I2 becomes 

ir = (^T-^A-^ep^e" • (64) 

We have proven that the DTT satisfies the Second Law for arbitrary values of the nonequi- 
librium tensor ^'^'^ , provided I2 is given by Eq. fl64l) . The entropy production is simply (recall 
Eq. iHD) 

St;, = D^'^^'^CpM, = ^^'"^P^ ■ (65) 
D. Exact hydrodynamic equations 

In this subsection we will obtain the explicit form of the equations of motion of the DTT. 
We first turn our attention to the quadratic part of the stress-energy tensor, T2. From 
Eq. fj44l) we get after some algebra 



(66) 



where ci = 2ci + C2. We note that [3^T2^ = [3yT2^ = and, of course (we calculated the 
coefficients for this to happen) rl*^ = 0. The divergence of T2'^ reads, 



3 (67) 
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From Eqs. (fTSi) . (133|) and (1671) . the conservation of the complete stress-energy tensor 
becomes 



' -1/ 



3 



+ (^^T^ - h.T-^C'^a^'^ (/3^A-- + p-'A^np^,, (68) 



where we used that p = aT'^ and p = p/3. 

In the spirit of divergence-type theories, the stress-energy tensor conservation should be 
supplemented with A^^'^ = J"'^, which stands on the same footing as the conservation equa- 



tions 



0,1 



38!]. We have already obtained 1°''^ in section UlI CI Together, they completely 
describe the space-time evolution of the system (within the hydrodynamic approximation). 

In this Section, we completed our first task of finding the generating function x and the 
source tensor J^''' that describes a conformal fiuid in fiat space-time. We constructed the 
DTT by requiring: (i) x is quadratic in deviations from equilibrium, represented by the 
dissipative tensor ^; (ii) the stress-energy tensor derived from x is traceless and has the 
correct conformal weight; (iii) the theory reproduces, at first order, the relativistic Navier- 
Stokes stress-energy tensor; and (iv) the theory satisfies the Second Law for arbitrary ^. We 
have also obtained the equations of motion of the exact theory, which will be used in Section 
IVII in the context of Bjorken expansion. 

IV. CAUSALITY 

In this section we investigate the causality properties of the DTT constructed above. As 
noted in Section [Tll causality is determined solely by the generating function x of the theory, 
and not by the source tensor 7°'^. In order to analyze causality of a DTT, let us define 

^A^^^M^M-^ and 

(69) 

M^j, = d'x/df3,dC^dC'' , 
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being {A, B) collective indices. The DTT is causal (in a set of fluid states near equilibrium) 
if the vector q5C,^5C,'~^ is time- 
from an equilibrium state 3^, 35, 36, 37, 38, 39, 40]. 



ike and future oriented for any displacements (5C ? "^C ) 



Since we are interested in proving causality for fluid states near equilibrium, it will be 
sufficient to deal with 



where 



,u 



E ^uS 



Ce ~ irri 5 ^ 

-LE 



0) 



denotes equilibrium values. 

The only non- vanishing terms of ^ 



Mi", 



0) 



are 



1) 



u5 



ML 



dl3^d(3-d^ 



and 



^E 



^E 



where X(o)y and X(2) are given by Eqs. ffT^ . fl7r|) and fIFIl) . respectively. 
Performing the corresponding derivatives, we get 



8an 



ML 



b 1^6^ + 65 - -g'^S^sj and 



where, for brevity, we have defined the tensor 



(70) 



(71) 



(72) 



(73) 



(74) 



Since Eqs. f[73|) are covariant, we can use any frame to study causality of our DTT. The 
frame = (1,0) turns out to be very convenient. In this frame we have 



SPu = {-t,w) 



(75) 
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(76) 



and 

^ A B^ 
B^ di 
B^ d2 

In writing the above equations, we have used the general decompositions of a vector and 
a tensor W'^'^ in time- and space-hke parts 



(77) 



with 



(78) 



V^Vu^ + X^" and 
W^"" = AuV + BV + u^C + E^'' 

V = -u^y^" , x^" - A/^v 

B^" = -A^W'^Ua 
C ^ -A/^W^Ua and 

Note that X^^u^ = B^u^ = C^u^ = E^^'^u^ = E^^Ui, = 0, and that, being real and sym- 
metric, E'^^ can be diagonahzed. In the case of S^^^" we have put 5^^^ — diag(cii, (i2, cia). 

— * 

Since 6^^''' should remain traceless and symmetric (as ^^^), we have B^ — — (0, B) and 
di + d2 + d3^ A. 

In the frame — (1)0) we have 



Mi's: 



SaT^ 



(79) 
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and 



2ci6q ( - 36,yi6TriSsjSej - 36,yi6eiSsjSTrj 



+ (c2 - ci)5^{5ui5so + 5uo55i)5ne 

Note that M^^^ is frame invariant (see Eq. ( !73l) ). 
Putting 



= mi:. 



= Mi:. 



dC^dp"" and 



we get, from Eqs. ([75]), ([76]), ([79]) and (IHOD . 



and 



= 8aT| Tst^ + ^ w,^ -2tw'j , 
^ i=i ^ 



= hSMC" = b\ -tA + B-w, -tB^ + djw^ 
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si" = ( ci[12G + —A^], -[ci{3dj +A)+ Ac2dj]B^ 



(80) 



(81) 



(82) 



(83) 



(84) 



where G = J^Ad-i)"^, j = (1,2,3) and no sum is imphed in the spatial part. So, we must 
now see whether y^* = + z'^ + is time-hke and future oriented, for arbitrary values of 
the fluctuations. It is clear that r'^ poses no problem; it is time-like and future oriented if 
a > 0. This was expected since corresponds to a perfect fluid. 
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Before analyzing the vectors z'^ and s^, it will be convenient to determine the constraints 
that the transversality of ^'^'^ imposes on them. We have found before that, in order for 
the DTT to be consistent, the nonequilibrium tensor had to be transverse. From a physical 
point of view, this meant that the heat flow and the bulk viscosity remain zero when the 
conformal fluid departs from equilibrium. This is a sensible result since we want the theory 
to remain conformal even in the presence of dissipation. Therefore, the condition we must 
impose is 

wie + sen = (Me + s(3,)sr = o (ss) 

where we have used that ^'^''l^; = 0. Using Eqs. (|75i) and (175]1 we get 

pf.Se" = (^-iA + w ■ B, -IB^ + djW^^ (86) 

with t = {t + Tj^^), being Te the equihbrium temperature of the fluid. Requiring 5^'' to be 
transverse we obtain 

tA = w ■ B and 

(87) 

tB^ = djW^ (no sum in j) 



or 

Using these results obtained from the transversality of e'^, we can rewrite the vector z'^ 
in a simple way: 

." = ^(^b). (89) 

Therefore, putting l'^ = z^ + we have 

f = ^A+^{12G+^-^A^) and 

^ (90) 



In order to prove that (Z°)^ > (/*)^ (i.e. that is time-like) it is convenient to reexpress l'^ 
in tensorial notation. From Eq. (!87l) we get 

B = -{d ■ w) and 

\ (91) 
Tr(d) = j-{w-d-w) , 
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where d = dij is the spatial part of S^f^,y, which we do not assume to be diagonal. Using 

these relations, becomes 

70 & ^ , ^ 26ci , ^ , ^,0 12ci , , , 
/ = w ■ a - w -\ — — -iw ■ a. ■ w) + „^ d : d and 



-'-E 



tT, 



d- w — J^^^ {w ■ d ■ w){d ■ w) 



(92) 



(3ci + 4c2) , ,2 



in 



{d' ■ w) 



where d : d stands for Tr([d]^). It is clear that if w = 0, l'^ is trivially time-like and future- 
oriented, provided Ci > 0, and therefore the theory is causal in this case. By continuity, 
the DTT will remain causal provided w is not too large. One can actually quantify this (at 
lowest order) by keeping linear terms in w and requiring that (/°)^ > I ■ I, but the resulting 
expression is not too illuminating. 

V. ADIABATIC EXPANSION IN VELOCITY GRADIENTS 

In this section, we set up a consistent adiabatic expansion of the DTT to compare with 



2Q|, 



2l|. As 



previous approaches based on conformal invariants, first put forward in Refs. 
stated in the Introduction, we will limit ourselves to Minkowski space-time. 

For a conformal fiuid in fiat space-time, the dissipative part of the stress-energy tensor 
complete at second-order in velocity gradients can be written as {20, [2I, 



30 



1 



'per 



-O 



3 



(93) 



where the subscript ci. is a remainder that this form of T2 is constructed from conformal 
invariants (as explained in detail in Refs. {20 . 



311]). D = u^d'^ is the convective 
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time derivative, (r^, Aj) are second-order transport coefficients, and VLp\ is the vorticity. As 
already mentioned, this expression for T2 represents an extension of Israel- Stewart entropy- 
wise approach. 

We will now show that, for the case A2,3 = 0, r^^ can be obtained from a consistent 
adiabatic expansion (at second-order in velocity gradients) of the exact hydrodynamic equa- 
tions. We start by requiring that calculated from x be equal to r^^^ calculated from 



22 



second-order (in velocity gradients) conformal invariants. We have 



+ — '-'(1) <7piipA + AsDj^j^j ^t'p^t'pX 



where we have put (see Eq. 



) 



Putting 



^(1) -(1) -(1) 

3 ^« 



r =eri)+e2) = -v^+er2) (96) 



in Eq. fl94l) . and retaining terms up to second order we get 



+ + (97) 



3-^(1) ^Ip^lpX - -^""^ ' 



Using Eq. (JHSD, we can rewrite the last equation more explicitly 

where we introduced < . . . > to denote the spatial, symmetric and traceless projection of a 
tensor: 

5<M-> = . (99) 

At second-order in velocity gradients, the equation A^.^"' = reads 

V 5/5. 5/5. 4p.J/^.;. + G ^^^^4.,;. ^^^^^ 
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where we have used Eqs. ( 155|) and (16^ . Exphcitly, we have 



r^a7p(Tf(2) _ ^a-^ _ '-^'^ c{l) n 



6 



The crucial point is that, in order for r''^ to be derivable from the DTT, both expressions 
for ^^^y given in Eqs. ( 1971) and fllOll) . should coincide. We see immediately that (actually, 
this equation holds for the exact ^pa) 



Using Eq. fl56l) . we can rewrite Eq. fllOip as 



+ 2(2ci + C2)T-' (^eS/^r. + eS/?^,) (103) 

The third term becomes 

- 12cir"^5f;]''"De^^^ = <D(t"^> , (104) 

which reproduces the first term of Eq. (1981) if = 12ci/(6T^). Using that 

u^g = -3D\nT , (105) 

it can be seen that the first term of Eq. (11031) reproduces the second term of Eq. (!98|) . The 
second term of Eq. fll03p reproduces the third and fourth terms of Eq. fl98l) . provided 

3 r 

Ai = ^ -4 + ^ (2C1 + C2) . (106) 

The last term of (I103P reproduces the fifth term of Eq. fl98l) if 2ci + C2 = ?)gT~'^rf. However, 
it is not possible to reproduce, from Eq. (I103p . the vorticity terms of Eq. (p8|) . So, we 
conclude that the DTT we have constructed is limited to the case A2.3 = 0. We note that 
this is not a serious restriction on the application of the DTT to heavy-ion collisions (see 
especially Ref. j^). 

We have proven that (for A2,3 = 0) r^^^, as given by Eq. ( I93l) . can be obtained from a 
consistent adiabatic expansion (at second-order in velocity gradients) of the exact divergence- 
type theory we have developed. This is one of the most important results of this work. 
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We have already proven that the DTT satisfies the Second Law. It is clear that its 
adiabatic expansion satisfies it too. It is interesting to remark that, when expanding the 
entropy production given in Eq. (1^ . i.e. when putting ^q,^ = ^a-y +^07, terms up to fourth 



order in velocity gradients arise. This agrees with the entropy production 
by Loganayagam in Ref. [30], based on the developments of Refs. [20, 



brm calculated 



2l| . Dropping 



fourth-order terms in the entropy production, although it may be justified under some 
circumstances, actually spoils the consistency of the adiabatic expansion (see Ref. 3l|] for 
interesting discussions on higher order terms in the entropy production). 



VI. BOOST INVARIANT FLOW 



We will now o 
(see also Refs 



btain the equations of motion of the DTT for the case of Bjorken flow 43 1 
, [13, [2^), which, besides of being much more simple than general flow. 



is a successful toy model of heavy-ion collisions in the mid-rapidity region. The comparison 
between the equations of the DTT and the second-order ones for the case of boost invariant 
flow is interesting because it clearly shows the difference between the exact and truncated 
equations, in a relatively simple situation. In the last part of this section, we compare the 
numerical solution to the exact and truncated equations. 

The motion in the Bjorken flow is a ID expansion, along an axis which we choose to be 
z, with local velocity equal to z/t. It is convenient to choose comoving coordinates (Milne 
coordinates), proper time r and rapidity ip, given by 

T = Vt"^ — z'^ and ip = arctanh(2;/t) . (107) 

The advantage of using these coordinates is that each element is at rest: (m"^, u^, u^) = 
(1,0,0). Although the velocity vector is constant, the dynamics is nontrivial because not 
every Christoffel symbol is zero. The metric tensor is 

fif^j, = diag(5(^^, g^^, gyy, g^^) = diag(l, -1, -1, -r^) (108) 

where {x, y) denote transverse directions, so we have 

D = u^d^ dr and 

1 (109) 

d,u^^ ^ - . 
r 
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The only nonvanishing component of the dissipative part of the stress-energy tensor is the 
{ip, ip) component. Note that the motion is irrotational, and that the energy density and the 
dissipative part of the stress-energy tensor only depend on proper time (i.e. are independent 
of rapidity). The only conservation equation that is nontrivial for Bjorken flow is the energy 
equation, i.e. m^T.^*^ = 0, where T'^'^ is the complete stress-energy tensor. 



^ , . . (110) 

± 



A. Second-order theory 

At second order in velocity gradients, the hydrodynamic equations for Bjorken flow are 
(see Refs. jsjls], Q for detailed discussions) 

d^p = + 5? with 

T T 

where, in the notation used here, 

nA'- = rr + rr. (Ill) 

Actually, the differential equation for showed in Eq. ( lllOp is exact up to terms which 
are second-order in velocity gradients. It is obtained by replacing a^'^ by U'^'^ in the gradient 
expansion of the latter [2(|] (see Eq. (l93l)). 

For a conformal perfect fluid in d=4, p(r) = Cr~^/^, where C is a constant. Due to 
conformal invariance, the viscosity and the second-order transport coefficients must scale as 
follows: 

/ \ 3/4 / N -1/4 / N 1/2 

where rjo, r° and A° are constants. 

Note that the Navier-Stokes equations are recovered formally by setting r^r, Ai — > 0, 
whereby 
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B. Divergence-type theory 



Projection of Eq. (168|) onto leads to 

Dp = -(p + p + h,T-'C^^a,] V.u'^ 

For Bjorken flow, energy conservation reads 

drP = -^(^P + P + ^(2ci + c,)T-Xr^ + kt (115) 
while the equation A^^'^ = 1°'^ becomes 

~ — ~ 3T7 ^ 73 



(116) 



where we made use of Eq. (11051) . 

The DTT as well as the second-order theory reduce to Eckart's theory when retaining 
first-order velocity gradients. Therefore, it is clear that the hydrodynamic equations of both 
theories must coincide in that limit (of course, this statement is valid for general flow, but 
we will discuss Bjorken flow only). Noticing that at first order in gradients we can write 
n^"^ = h^l^^ and (ci, C2, (7, Ttt, Ai) = 0, we immediately see that the hydrodynamic equations 
of the DTT and the second-order theory (Eqs. (I115|116l) and Eqs. (II 101) . respectively) 
coincide. From Eq. (11161) we recover the Navier-Stokes limit given in Eq. (11131) . 

The comparison between the equations of both theories beyond first order in velocity 
gradients becomes quite complicated, because, being r|"^ quadratic in , 11^" and are 
not linearly related anymore (see Eq. (Illip ). In order to carry out the this comparison, 
we solve both sets of differential equations numerically in the next subsection. Before doing 
that, it is convenient to reexpress (using the results of the previous section) the equations 
of the DTT in terms of (77, r^, Ai) instead of (6, Ci, C2). Without loss of generality, we can fix 
h = 7] (this means ^^^^ = — cr'^''), whereby 

— — T and Ci = 

3 V 2 r/2 y 12 



C2 = -V(^ + ^) and ci = ^. (117) 



Note that g is completely specified once Ci and C2 are known: 

_ AiT^ 

^~ V ■ 
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(118) 



Equations (11151) and (11161) then read 



with 



and 



with 



drP = -Up + p-TF,[^^iA+y^ (119) 



drd = + E,d + (121) 





At 


E2 = 


-(^ 


Es = 







r2 



(122) 



The dissipative part of the stress-energy tensor in the DTT is constructed from the solution 
to Eq. ffml) . We have 

+ < = < + i^iT[eJf . (123) 



C. Comparison of numerical solutions 

In this section, we compare the solutions to the hydrodynamic equations of the DTT, 
second-order and Navier-Stokes theories. We focus on the inverse Reynold's number 

R-' = ^, (124) 
P + P 

and on the pressure isotropy 

^ - ^ . (125) 

Pt p + n*/2 

These two quantities are relevant parameters to characterize the hydrodynamic evolution 



fsee, for instance, Refs. 



an 



13 



14 



ISip. Ideal fluids are characterized by i? ^ = and 



Pl/Pt = 1- Note that, as already mentioned, in the DTT we have 11^ = rf^ + r^. When 
solving the hydrodynamic equations, one must bear in mind that the transport coefficients 
are functions of the energy density p, as given by Eq. (I112p . In particular, we will focus on 
the strongly-coupled SYM plasma, for which we have [20] 

r^ = 2(2-ln2)^ and Ai = ^ , (126) 
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FIG. 1: Inverse Reynold's number i? as a function of proper time, for the DTT, second-order 
and Navier-Stokes theories with rj/s = 0.09. 

where s is the entropy density. 

In the following, we present the results for two relevant values of r]/s. This value is 
modified by changing the value of r/o- We consider rj/s = 0.09, which is very close to 
the lower bound imposed by the AdS/CFT correspondence {t]/s > 1/Att), and rj/s = 0.375, 
which is close to the upper bound for the quark-gluon plasma found by comparing dissipative 
hydrodynamics to elliptic flow measurements {rj/s < 0.5) [8]. As initial conditions, we set 
n^(ro) = 0, To = 0.5 fm/c and p{tq) = 10 GeV/fm^ in all calculations. 

In Figured] we compare the evolution of the inverse Reynold's number with proper time 
for the DTT, the second-order and Navier-Stokes theories with rj/s = 0.09. The most 
important feature is that the DTT shows a faster approach to ideal hydrodynamics. Figure 
|2] shows the same comparison but for i]/s = 0.375. As in the previous case, the DTT shows 
a faster approach to the ideal fluid behaviour. 

In Figure [3] we show the evolution of the pressure isotropy for rj/s = 0.09. It is clearly 
seen that the approach to ideal hydrodynamics is faster in the DTT, which also occurs with 
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FIG. 2: Inverse Reynold's number i? ^ as a function of proper time, for the DTT, second-order 
and Navier-Stokes theories with r]/ s = 0.375. 

r]/s = 0.375 (Figure SD. 

We note that, with respect to the second-order theory, our results are in good agreement 



with those of previous studies 



ll|,[l2 



13 



14 



15|. Taking into account the behaviour of the 



two quantities that we analyzed, we arrive at the important conclusion that the relaxation 
towards ideal hydrodynamics is faster in the DTT than in the second-order theory. This 
means that, as expected on theoretical grounds, the hydrodynamic evolution in the DTT is 
closer to that obtained from transport theory (see in particular the detailed comparison be- 
tween Navier-Stokes, Israel- Stewart and covariant transport theory carried out by Houvinen 
and Molnar in Ref. [l4|). 



VII. SUMMARY AND CONCLUSIONS 



In this work, we have studied the (nonlinear) hydrodynamical description of a conformal 
field within the theoretical framework of divergence-type theories. We proved that the theory 
we develop is causal (in a set of fluid states near equilibrium) and satisfies the Second Law 
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FIG. 3: Pressure isotropy Pl/Pt as a function of proper time, for the DTT, second-order and 
Navier-Stokes theories with r]/s = 0.09. 

exactly. Since it does not rely on gradient expansions, it goes beyond second-order (in 
velocity gradients) theories, thus being a closed theory. However, it is limited to the case 
where the second-order transport coefficients A2 and A3 vanish. For this case, we showed 
that the second-order stress-energy tensor constructed from conformal invariants 0,0, Q 
can be consistently derived via an adiabatic expansion from the DTT. 

As the most simple example, we have also obtained the hydrodynamic equations of the 
DTT for Bjorken flow, and compared them, analytically and numerically, with those of 
second-order and Navier-Stokes theories. The numerical calculations indicate that the re- 
laxation towards ideal hydrodyanamics is substancially faster in the DTT as compared to 
the second-order theory. This indicates that the DTT is a better approximation to transport 
theory than the second-order theory, as expected since the former includes all-order velocity 
gradients. 

As stated in the Introduction, we think that the theory we have presented may be useful 
in the analysis of early-time dynamics and in the evolution of initial state fluctuations in 
heavy-ion collisions, essentially because the theory is not based on an expansion in velocity 
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FIG. 4: Pressure isotropy Pl/Pt as a function of proper time, for the DTT, second-order and 
Navier-Stokes theories with r]/s = 0.375. 

gradients. The extension of the DTT to include the case A2,3 7^ is also interesting. Work 
is in progress along these lines. 
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